arXiv:1502.00794vl [cs.IT] 3 Feb 2015 


EURASIP JOURNAL ON WIRELESS COMMUNICATIONS AND NETWORKING, SPECIAL ISSUE ON FEMTOCELL NETWORKS, 2010 


1 


Best Signal Quality in Cellular Networks: 
Asymptotic Properties and Applications to Mobility 
Management in Small Cell Networks 

Van Minh Nguyen, Francois Baccelli, Laurent Thomas, and Chung Shue Chen 


Abstract —The quickly increasing data traffic and the user 
demand for a full coverage of mobile services anywhere and 
anytime are leading mobile networking into a future of small 
cell networks. However, due to the high-density and randomness 
of small cell networks, there are several technical challenges. 
In this paper, we investigate two critical issues: best signal 
quality and mobility management. Under the assumptions that 
base stations are uniformly distributed in a ring shaped region 
and that shadowings are lognormal, independent and identically 
distributed, we prove that when the number of sites in the 
ring tends to infinity, then (i) the maximum signal strength 
received at the center of the ring tends in distribution to a 
Gumbel distribution when properly renormalized, and (ii) it 
is asymptotically independent of the interference. Using these 
properties, we derive the distribution of the best signal quality. 
Furthermore, an optimized random cell scanning scheme is 
proposed, based on the evaluation of the optimal number of 
sites to be scanned for maximizing the user data throughput. 

Index Terms —Small cell networks, maximum SINK, handover, 
random cell scanning, extreme value theory. 

I. Introduction 

Mobile cellular networks were initially designed for voice 
service. Nowadays, broadband multimedia services (e.g., video 
streaming) and data communications have been introduced into 
mobile wireless networks. These new applications have led to 
increasing traffic demand. To enhance network capacity and 
satisfy user demand of broadband services, it is known that 
reducing the cell size is one of the most effective approaches 
0-0 to improve the spatial reuse of radio resources. 

Besides, from the viewpoint of end users, full coverage 
is particularly desirable. Although today’s macro and micro 
cellular systems have provided high service coverage, 100%- 
coverage is not yet reached because operators often have many 
constraints when installing large base stations and antennas. 
This generally results in potential coverage holes and dead 
zones. 

A promising architecture to cope with this problem is that 
of small cell networks 0 , 0 . A small cell only needs 
lightweight antennas. It helps to replace bulky roof top base 
stations by small boxes set on building facade, on public 
furniture or indoor. Small cells can even be installed by end 
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users (e.g., femtocells). All these greatly enhance network 
capacity and facilitate network deployment. Pervasive small 
cell networks have a great potential. For example, Willcom 
has deployed small cell systems in Japan 0, and Vodafone 
has recently launched home 3G femtocell networks in the UK 
0 . 

In principle, high-density and randomness are the two basic 
characteristics of small cell networks. First, reducing cell size 
to increase the spatial reuse for supporting dense traffic will 
induce a large number of cells in the same geographical 
area. Secondly, end users can set up small cells by their 
own means ED- This makes small cell locations and coverage 
areas more random and unpredictable than traditional mobile 
cellular networks. The above characteristics have introduced 
technical challenges that require new studies beyond those 
for macro and micro cellular networks. The main issues 
concern spectrum sharing and interference mitigation, mobility 
management, capacity analysis, and network self-organization 
0, 0- Among these, the signal quality, e.g., in terms of 
signal-to-interference-plus-noise ratio (SINR), and mobility 
management are two critical issues. 

In this paper, we first conduct a detailed study on the 
properties of best signal quality in mobile cellular networks. 
Here, the best signal quality refers to the maximum SINR 
received from a number of sites. Connecting the mobile to 
the best base station is one of the key problems. The best 
base station here means the base station from which the 
mobile receives the maximum SINR. As the radio propagation 
experiences random phenomena such as fading and shadowing, 
the best signal quality is a random quantity. Investigating 
its stochastic properties is of primary importance for many 
studies such as capacity analysis, outage analysis, neighbor 
cell scanning, and base station association. However, to the 
best of our knowledge, there is no prior art in this area. 

In exploring the properties of best signal quality, we focus 
on cellular networks in which the propagation attenuation 
of the radio signal is due to the combination of a distance- 
dependent path-loss and of lognormal shadowing. Consider a 
ring B of radii li,mn and Rb such that 0 < J? m in < Rb < oo. 
The randomness of site locations is modeled by a uniform 
distribution of homogeneous density in B. Using extreme 
value theory (c.f., 0, 0), we prove that the maximum 
signal strength received at the center of B from n sites in 
B converges in distribution to a Gumbel distribution when 
properly renormalized and it is asymptotically independent of 
the total interference, as n —> oo. The distribution of the best 
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signal quality can thus be derived. 

The second part of this paper focuses on applying the 
above results to mobility support in dense small cell networks. 
Mobility support allows one to maintain service continuity 
even when users are moving around while keeping efficient 
use of radio resources. Today’s cellular network standards 
highlight mobile-assisted handover in which the mobile mea¬ 
sures the pilot signal quality of neighbor cells and reports the 
measurement result to the network. If the signal quality from 
a neighbor cell is better than that of the serving cell by a 
handover margin, the network will initiate a handover to that 
cell. The neighbor measurement by mobiles is called neighbor 
cell scanning. Following mobile cellular technologies, it is 
known that small cell networking will also use mobile-assisted 
handover for mobility management. 

To conduct cell scanning m-m, today’s cellular net¬ 
works use a neighbor cell list. This list contains information 
about the pilot signal of selected handover candidates and 
is sent to mobiles. The mobiles then only need to measure 
the pilot signal quality of sites included in the neighbor cell 
list of its serving cell. It is known that the neighbor cell 
list has a significant impact on the performance of mobility 
management, and this has been a concern for many years in 
practical operations lfl3ll . fl4l| as well as in scientific research 
m-m. Using neighbor cell list is not effective for the 
scanning in small cell networks due to the aforementioned 
characteristics of high-density and randomness. 

The present paper proposes an optimized random cell 
scanning for small cell networks. This random cell scanning 
will simplify the network configuration and operation by 
avoiding maintaining the conventional neighbor cell list while 
improving user’s quality-of-service (QoS). It will also be 
implementable in wideband technologies such as WiMAX and 
LTE. 

In the following, SectionHIldescribes the system model. Sec¬ 
tion |III] derives the asymptotic properties and the distribution 
of the best signal quality. Section IV presents the optimized 
random cell scanning and numerical results. Finally, Section V 
contains some concluding remarks. 

II. System Model 

The underlying network is composed of cells covered by 
base stations with omni-directional antennas. Each base station 
is also called a site. The set of sites is denoted by il C N. 
We now construct a model for studying the maximum signal 
strength, interference, and the best signal quality, after spec¬ 
ifying essential parameters of the radio propagation and the 
spatial distribution of sites in the network. 

As mentioned in the introduction, the location of a small 
cell site is often not exactly known even to the operator. The 
spatial distribution of sites seen by a mobile station will hence 
be treated as completely random fl9l and will be modeled by 
an homogeneous Poisson point process (20) with intensity A. 

In the following, it is assumed that the downlink pilot signal 
is sent at constant power at all sites. Let f? m ; n be some strictly 
positive real value. For any mobile user, it is assumed that the 
distance to his closest site is at least R m i n and hence the path 


loss is the far-held. So, the signal strength of a site i received 
by a mobile at a position y £ R 2 is given by 

Pi(y) = MW ~ X i|) _/3 2Q, for |y — x,| > i? min , (1) 


where x, <£ M 2 is the location of site i, A represents the 
base station’s transmission power and the characteristics of 
propagation, /3 is the path loss exponent (here, we consider 
2 < /3 < 4), and the random variables AQ = 10 x * dB / 10 , 
which represent the lognormal shadowing, are defined from 
{X dB , i = 1,2,...}, an independent and identically dis¬ 
tributed (i.i.d.) sequence of Gaussian random variables with 
zero mean and standard deviation <TdB- Typically, o,\\> is 
approximately 8 dB ED, E2- Here, we consider that fast 
fading is averaged out as it varies much faster than the 
handover decision process. 

Cells sharing a common frequency band interfere. Each cell 
is assumed allocated no more than one frequency band. Denote 
the set of all the cells sharing frequency band fc-th by ilk, 
where k = 1 So ilk fl ilk' = 0 for k ^ k' , and 

Uf = i Hk — il- The SINR received at y £ R 2 from site i £ ilk 
is expressible as 


C»(y) 


_ Pj(y) _ 

No + p j(y) ’ 


for i £ ilk, 


( 2 ) 


where Nq is the thermal noise average power which is assumed 
constant. For notational simplicity. Let A := A/Nq. Then 

Ci(y) is g iven by 


C»(y) = 


p i{ y) 


i + E^ije s-4 p j( y) ’ 

In the following, we will use (0 instead of 0. 


for i £ ilk • (3) 


III. Best Signal Quality 

In this section, we derive the distribution of the best signal 
quality. Given a set of sites S C il, the best signal quality 
received from S' at a position y £ R 2 , denoted by Y$(y), is 
defined as: 

F s (y) = maxQ(y). (4) 


Let us first consider a single-frequency network (i.e., K = 

1 ). 

Lemma 1: In the cell set S of single-frequency network , 
the site which provides a mobile the maximum signal strength 
will also provide this mobile the best signal quality, namely 


Ys{ y) = 


M s (y) 

1 + T(y)-Msr(y)’ 


VyeR 2 , 


(5) 


where 

Ms{y) = max Pi (y) 


is the maximum signal strength received at y from the cell set 
S, and 

f(y) = E p *(y) 

ie n 


is the total interference received at y. 

Proof: Since Q(y) = Pi( y)/{l + I(y) ~ Pi{ y)} and 
p(y) < /(y), 0 follows from the fact that no matter which 
cell i £ il is considered, J(y) is the same and from the fact 
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that x/(c — x ) with c constant is an increasing function of 
x < c. □ 

Let us now consider the case of multiple-frequency net¬ 
works. Under the assumption that adjacent-channel interfer¬ 
ence is negligible compared to co-channel interference, cells 
of different frequency bands do not interfere one another. Thus, 
for a given network topology T, the SINRs received from cells 
of different frequency bands are independent. In the context 
of a random distribution of sites, the SINRs received from 
cells of different frequency bands are therefore conditionally 
independent given T. Write cell set S as 

K 

s= LJ{Sfc : S k cn k }, 

k =1 

with Sk the subset of S allocated to frequency k. Let 


Y S k { y) = maxCi(y) 

lESk 


be the best signal quality received at y from sites which 
belong to S k - The random variables {Yg k (y),k = 1 
are conditionally independent given T. As a result, 

K 

P{^(y) < 7 1T} = I] (y) < 7 1 T}. 

k= 1 

Remark. For the coming discussions, we define 


I s(y) = J2 Pi ( y ) 

ies 


which is the interference from cells in set S. In the following, 
for notational simplicity, the location variable y appearing in 
Yg( y), Mg{ y), Ig{ y) and /(y) will be omitted in case of no 
ambiguity. We will simply write Ys, Mg, Is, and I. Note that 
I S <I since S C SI. 

Following Lemma [I] the distribution of Yg can be deter¬ 
mined by the joint distribution of Mg and I, which is given 
below. 

Corollary 1: The tail distribution of the best signal quality 
received from cell set S is given by 

roc riff-U-l 

f y s ( 7 ) = / / f(i,Ms) (v, u)dvdu (6) 

J u= 0 Jv=u 

where A ^ le joint probability density of I and Mg. 

Proof: By Lemma Q] we have 


P{*s > 7 } = P{M S /(1 +1 - M s ) > 7 } 
= P{/< —Ms-1} 


u=0 J v=u 


f{i,M s )( v y u )dvdu. 


□ 

In view of Corollary [Q we need to smdy the properties 
of the maximum signal strength Mg as well as the joint 
distribution of Mg and I. As described in the introduction, 
in dense small cell networks, there could be a large number 
of neighbor cells and a mobile may thus receive from many 
sites with strong enough signal strength. This justifies the use 
of extreme value theory within this context. 


For some /i’ m i„ and Up such that 0 < i? m i n < Rb < 00 , let 
6 C l 2 be a ring with inner and outer radii Ii ]n \ n and Rb, 
respectively. In this section, we will establish the following 
results: 

(i) The signal strength Pi received at the center of B 
belongs to the maximum domain of attraction (MDA) 
of the Gumbel distribution (c.f., Theorem |T] in Sec¬ 
tion QTLA]i. 

(ii) The maximum signal strength and the interference re¬ 
ceived at the center of B from n sites therein are 
asymptotically independent as n —> 00 (c.f., Corollary 0 
in Section UlI-AI) . 

(iii) The distribution of the best signal quality is derived (c.f.. 
Theorem |2] in Section Illl-Cb . 

A. Asymptotic Properties 

To begin with, some technical details need to be specified. 
Given a ring B as previously defined, we will study metrics 
(such as e.g., signal strength, interference, etc.) as seen at the 
center of B for a set S C Cl of n sites located in B. We will 
use the notation M n , Y n , and /,, instead of Mg, Yg, and Ig, 
respectively, with 

n 

M n = max Pi, I n = Y'' Pi, Y n = max Q. 

i=l,i£S ^ i=l, i€S 

2=1,z£«S 

Lemma 2: Assume that 0 < R m in < Rb < 00 , that sites 
are uniformly distributed in B, and that the shadowing Xi 
follows a lognormal distribution of parameters (0, ax)- Then 
the cdf of the signal strength Pi received at the center of B 
from a site located in B is given by: 

Fp{x) = c{a~^ Gi(x) — b~^G 2 {x) 

— e v x~^G^,{x) + e l 'x~^G4(x)J (7) 

where a^= AR~ P , b = AR~f n , c = A? (i?| - in )~\ 

v = 2a\//3 2 , and Gj, j = refers to the cdf of a 

lognormal distribution of parameters ( Pj,crx), in which 

pi = log a, p 3 = pi + 2 ax/ft, 

P 2 = log b, PA = P 2 + 2.a 2 x /p. 

Proof: See Appendix lAl □ 

Under the studied system model, {Pi, i = 1,2,...} are 
independent and identically distributed (i.i.d.), and so the cdf 
Fm„ and probability density function (pdf) Jm u of M n are 
directly obtained as follows: 

Corollary 2: Under the conditions of Lemma [2] the cdf and 
the pdf of M n are given respectively by: 

F Mn {x) =*p(x), ( 8 ) 

/m„ (x) = nf P (x)Fp~ 1 (x), (9) 

where Fp{x) is given by ( 0 , and fp is the pdf of Pi, fp{x) = 
dFp(x)/dx. □ 

Since M n is the maximum of i.i.d. random variables, we can 

also study its asymptotic properties by extreme value theory. 
Fisher and Tippett [9l Thm. 3.2.3] proved that under appropri¬ 
ate normalization, if the normalized maximum of i.i.d. random 
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variables tends in distribution to a non-degenerate distribution 
H, then H must have one of the three known forms: Frechet, 
Weibull, or Gumbel distribution. In the following, we prove 
that Pi belongs to the MDA of a Gumbel distribution. First of 
all, we establish the following result that is required to identify 
the limiting distribution of M n . 

Lemma 3: Under the conditions of Lemma [2] the signal 
strength received at the center of B from a site located in B 
has the following tail equivalent distribution: 


exp ( — (log x - /i 2 ) 2 /(2 a 2 x )) 
PX ^ (log x — p 2 ) 2 / (2<7jt) 
2y/2TTo x G 2 (x) 

~ k —^-, as x —> oo, 

log X- H 2 


(10a) 

(10b) 


where G 2 {x) = 1 - G 2 (x), and k = 

V2ir£ H s -Hmin _ 

Proof: See Appendix IB I □ 

Equation (11 Obb shows that the tail distribution of the signal 
strength Pi is close to that of G 2 , although it decreases more 
rapidly. The fact that G 2 determines the tail behavior of Fp 
is in fact reasonable, since G 2 is the distribution of the signal 
strength received from the closest possible neighboring site 
(with b = AR~? n and ax)- The main result is given below. 

Theorem 1: Assume that 0 < R m < Rp < oo, that sites 
are uniformly distributed in B, and that shadowings are i.i.d. 
and follow a lognormal distribution of parameters (0,<7x) 
with 0 < ax < oo. Then there exists constants c n > 0 and 
d n £ R such that: 


c„ 1 (M. n — d n ) —> A as n —> oo, 


( 11 ) 


where A is the standard Gumbel distribution: 


A(x) = exp{—e x }, 

and 7 represents the convergence in distribution. A possible 
choice of c n and d n is: 


Cn 

dn 


cr x (21ogn) id n , 

exp {p -2 + o x (\/2 log n + 


— log log n + log k 
y/2 logn 


)}. 

( 12 ) 


with k given by Lemma [3 

Proof: See Appendix© □ 

By Theorem Q] the signal strength belongs to the MDA 
of the Gumbel distribution, denoted by MDA(A). From |[23l . 
[24J, we have the following corollary of Theorem |T| 

Corollary 3: Let a p be the variance and pp be the mean 
of signal strength Pp Let I n = ( I n — npp)/(y/nop). Let 
M n = (M n — d n )/c n , where c n and d n are given by (112b . 
Under the conditions of Theorem [7] 


(M n ,I n )A (A, $) as n -»■ oo, (13) 


where A is the Gumbel distribution and $ the standard Gaus¬ 
sian distribution, and where the coordinates are independent. 

Proof: Conditions 0 < R min and ox < oo provide a 2 P < 
vai'{ AR~Xj } < oo. Then the result follows by Theorem Q] 

and (231, EH- □ 


Note that the total interference / can be written as I = 
In + In where denotes the complement of I n in /. Under 
the assumptions that the locations of sites are independent 
and that shadowings are also independent, I n and are 
independent. The asymptotic independence between M n and 
I n thus induces the asymptotic independence between M n and 
I. This observation is stated in the following corollary. 

Corollary 4: Under the conditions of Theorem [7] M n and 
I are asymptotically independent as n —> oo. □ 

This asymptotic independence facilitates a wide range of 
studies involving the total interference and the maximum 
signal strength. This result will be used in the coming sub¬ 
section to derive the distribution of the best signal quality. 

Remark'. The asymptotic properties given by Theorem Q] 
and Corollaries [3] and [4] hold when the number of sites in a 
bounded area tends to infinity. This corresponds to a network 
densification process in which more sites are deployed in 
a given geographical area in order to satisfy the need for 
capacity, which is precisely the small cell setting. 

B. Convergence Speed of Asymptotic Limits 

Theorem Q] and Corollaries [3 and Q] provide asymptotic 
properties when n —> oo. In practice, n is the number of cells 
to be scanned, and so it can only take moderate values. Thus, 
it is important to evaluate the speed of convergence speeds 
of (fill and ( 1 1 3b . We will do this based on simulations and 
will measure the discrepancy using a symmetrized version of 
the Kullback-Leibler divergence (the so called Jensen-Shannon 
divergence (JSdiv)). 

Let us start with some numerical evaluations of the con¬ 
vergence of M n to its limiting distribution. Fig. |T] a) shows 
Fy: ; for different n and compares to empirical simulation 
results. As expected the analytical distributions obtained by 
© of Corollary [3 match with the empirical distributions 
for all n. Fig. [TJb) plots the analytical distribution and its 
limiting distribution, i.e. the Gumbel distribution A. There is 
a discrepancy in the negative regime (see the circled region in 
FigdJb)). It is worth pointing out that as a maximum of signal 
strengths, M n > 0 and thus M n > —d n /c n = —yj2\ogn/ax 
since M n = ( M n — d n )/c n . This means that F^ In (x) = 0, 
Vx < — \J2 log n/ox, whereas A(x) > 0, Vir > — oo. This 
explains the gap observed for n small. This dissimilarity 
should have limited impact as long as we only deal with 
positive values of M n (respectively, M n > —\J 2 \ogn/ax)- 

We now study the symmetrized divergence between the 
analytical and limiting distributions of M n for some moderate 
values of n and under different Ode and /3. The convergence 
is best for OdB around 10 dB and /3 around two to four. For 
practical systems, OdB is approx. 8 dB and 2 < P < 4. 
We compute the Jensen-Shannon divergence for (3 = 3 and 
<7dB = 8 and plot the results in Fig. [3 a) and (b), respectively. 
For these (and other) values (within the range given above) of 
<7dB and if M n and A have low divergence. 

Let us now measure the (dis)similarity between the empiri¬ 
cal joint distribution, P {M n < u, l n < v}, and the product of 
the empirical marginal distributions, P {M n < u) x P{1„ < 
v}. Fig.|3shows an example with n = 50, /? = 3 and oas = 8. 
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normalized signal strength [decimal scale] 


o Emperical n=5 

.Analytical n=5 

A Emperical n= 10 
- Analytical n=10 

“ Emperical n=10 2 

2 

- Analytical n=10“ 


(a) Comparison between analytical and empirical CDFs 


Fig. 1. CDF of M„ under different n: ct,\ g = 8, /3 = 3. 



normalized signal strength [decimal scale] 

(b) Comparison between analytical and limiting CDFs 



We see that these two density functions are very similar. Fig. [4] 
compares these two density functions for different values of 
OdB and (3. Within the range defined above, the divergence 
between the two distributions is again small. Comparing Fig. [2] 
and Fig. [4] one can conclude that even if the convergence 
of M r , 4 A remains slow, M n and I n quickly become 
uncorrelated. Thus, the independence between M n and /,, 
holds for moderate values of n , i.e., 

f(M n ,i n )(u,v) w f Mn (u ) x fi n (v), (14) 

and so the same conclusion holds for the independence be¬ 
tween M n and I. 


C. Distribution of the Best Signal Quality 

From the above results, we have the distribution of M n and 
the asymptotic independence between M n and I. In order to 
derive the distribution of the best signal quality, we also need 
the distribution of the total interference. 


Lemma 4: Assume that shadowings are i.i.d. and follow 
a lognormal distribution of parameters (0 ,<Jx), (3 > 2, and 
that sites are distributed according to a Poisson point process 
with intensity A outside the disk of radius f?. m i n > 0. Let I 
be the interference received at the disk center, and 4>i be the 
characteristic function of I. Then: 

1 ) 


<j)i{w) = exp{ — TtAa(A\w\) a 


A m 



1 - fx (sign(w)f) 
t a+l 


d t}, 
(15) 


where a = 2 /(3, and <px is the characteristic function 

ofXi. 

2) \(fi{w)\ p £ L for all p = 1,2,..., where L is the space 
of absolutely integrable functions. 

3) If is large, then (j)j admits the following approx¬ 

imation: 


(f I {w) ~ exp ( — <5|ui|“[l — jsign(tu) tan(—)]), (16) 
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f V, Xf T 

norm. M norm. I 


(norm. M , norm. I ) 



Fig. 3. Example of joint densities of M n and I n : n = 50, <t<jb = 8, /3 = 3. Here, norm. M n refers to M„, while norm. I n refers to I n . 




(a) Under different crj b: ft = 3 (b) Under different ft : tr^B = 8 

Fig. 4. Jensen-Shannon divergence between j , and fyy. X fj^. 


where 8 = 7tAA“ exp(ia 2 cr^.)r(l — a) cos(7tq!/2), with 
r(-) denoting the gamma function. 

Proof: See Appendix iDl □ 

Theorem 2: Under the assumptions of Lemma [7] let B be 
the ring of inner and outer radii R m i n and Rb, respectively. 
Denote the best signal quality received at the center of B from 
n sites in B by Y n . Assume that 0 < ax < oo, 0 < R min < 
Rb < oo, that AR~f m is large, and that ttX(R% — -R^in) > 
with high probability, where n is some positive integer. Then 
the tail distribution of Y n can be approximated by: 

r r ( ~*~ i o 

Fy n ( 7 ) ~ / { k {u) - e~ Su,a sin ( w —2) 

J-y L J o ™ 27 

x cos (wu + w - — — Sw a tan —-)dw;Idw. 

v 27 2 ' / 

Proof: See Appendix [E] 


The approximation proposed in Theorem 2 will be used in 
Section [IV] below. It will be validated by simulation in the 
context considered there. 

IV. Random Cell Scanning for Data Services 

In this section, the theoretical results developed in Sec¬ 
tion |III] are applied to random cell scanning. 

A. Random Cell Scanning 

Wideband technologies such as WiMAX, WCDMA, and 
LTE use a predefined set of codes for the identification of 
cells at the air interface. For example, 114 pseudo-noise 
sequences are used in WiMAX l25l , while 504 physical cell 
(17) identifiers are used in LTE |26j . When the mobile knows the 
identification code of a cell, it can synchronize over the air 
□ interface and then measure the pilot signal quality of the cell. 
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Therefore, by using a predefined set of codes, these wideband 
technologies can have more autonomous cell measurement 
conducted by the mobile. In this paper, this identification code 
is referred to as cell synchronization identifier (CSID). 

In a dense small cell network where a large number of cells 
are deployed in the same geographical area, the mobile can 
scan any cell as long as the set of CSIDs used in the network 
is provided. This capability motivates us to propose random 
cell scanning which is easy to implement and has only very 
few operation requirements. The scheme is detailed below: 

(1) When a mobile gets admitted to the network, its (first) 
serving cell provides him the whole set of CSIDs used 
in the network. The mobile then keeps this information 
in its memory. 

(2) To find a handover target, the mobile randomly selects a 
set of n CSIDs from its memory and conducts the stan¬ 
dardized scanning procedure of the underlying cellular 
technology, e.g., scanning specified in IEEE 802.16 H25I . 
or neighbor measurement procedure specified in 3G l l27l 
and LTE ifTZl . 

(3) The mobile finally selects the cell with the best received 
signal quality as the handover target. 

In the following, we determine the number of cells to be 
scanned which maximizes the data throughput. 

B. Problem Formulation 

The optimization problem has to take into account the two 
contrary effects due to the number of cells to be scanned. 
On one hand, the larger the set of scanned cells, the better 
the signal quality of the chosen site, and hence the larger the 
data throughput obtained by the mobile. On the other hand, 
scanning can have a linear cost in the number of scanned cells, 
which is detrimental to the throughput obtained by the mobile. 

Let us quantify this using the tools of the previous sections. 

Let W be the average cell bandwidth available per mobile 
and assume that it is a constant. Under the assumption of 
additive white Gaussian noise, the maximum capacity £„ that 
the mobile can have by selecting the best among n randomly 
scanned cells is 


4 = wao g (i + y„). (is) 

Hence 

E(U = IUE{log(l + Y n )} 

pOO 

= W log(l +7)/w(7)d7> 

J 7=0 

where fy n is the pdf of Y n . By an integration by parts of 
log(l + 7) and fy n ( 7)d7 = — dFy n (7), this becomes: 

Fy h) 

EU 4 = W / -j^d 7 . ( 19 ) 

4=0 1 + 7 

Note that E{£„} is the expected throughput from the best 
cell. Since Y n is the maximum signal quality of the n cells, 
Y n increases with n and so does £„■ Hence, the mobile should 
scan as many cells as possible. However, on the other hand, 
if scanning many cells, the mobile will consume much time 
in scanning and thus have less time for data transmission 


with the serving cell. A typical situation is that where the 
scanning time increases proportionally with the number of 
cells scanned and where the data transmission is suspended. 
This for instance happens if the underlying cellular technology 
uses a compressed mode scanning, like e.g., in IEEE 802.16 
l25l and also inter-frequency cell measurements defined by 
3GPP H2, l27l . 

Another scenario is that of parallel scanning-transmission: 
here scanning can be performed in parallel to data transmission 
so that no transmission gap occurs; this is the case in e.g., 
intra-frequency cell measurements in WCDMA | |27] and LTE 

D3- 

Let T be the average time during which the mobile stays 
in the tagged cell and receives data from it. Let s be the time 
needed to scan one cell (e.g., in WCDMA, the mobile needs 
s = 25 ms if the cell is in the neighbor cell list and s = 800 ms 
if not ll28l . whereas in WiMAX, s = 10 ms, i.e., two 5-ms 
frames). Let L(n) be the duration of the suspension of data 
transmission due to the scanning of the n cells: 


L{n) 


s x n if compressed mode is used, 

0 if parallel scan.-trans. is enabled. 


( 20 ) 


Finally, let E{£o} be the average throughput received from the 
serving cell when no scanning at all is performed (this would 
be the case if the mobile would pick as serving site one of the 
sites of set S at random). 

The gain of scanning n cells can be quantified by the 
following metric, that we will call the acceleration: 


T- E{U 

T ‘ E{£o} + L(n) ■ E{£ 0 } 
T E{U 
T + L(n) E{£ 0 } 


( 21 ) 


In this definition, T ■ E{£ n } (resp. T ■ E{£o} + L(n ) ■ E{£o})) 
is the expected amount of data transmitted when scanning n 
cells (resp. doing no scanning at all). We aim at finding the 
value of n that maximizes the acceleration p n . 

It is clear that T/(T + L(n)) = 1 when (i) T — > 00 , i.e., 
the mobile stays in and receives data from the tagged cell 
forever, or (ii) L(n) = 0, i.e., parallel scanning-transmission 
is enabled. In these cases, p n increases with n and the mobile 
“should” scan as many cells as possible. However, p n is 
often concave and the reward of scanning then decreases. To 
characterize this, we introduce a growth factor g defined as 
follows: 

a_ Pn _ T + L(n — 1) E{4} 

9n ~ Pn-1~ T + L(n) ’ 


Special cases as those considered above can be cast within 
a general framework which consists in finding the value of 
n that maximizes p n under the constraint that g n > 1 + A g , 
where A 9 > 0 is a threshold. 


[n this mode, scanning intervals, where the mobile temporarily suspends 
data transmission for scanning neighbor cells, are interleaved with intervals 
where data transmission with the serving cell is resumed. 
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Fig. 5. Numerical results in the random cell scanning optimization. 
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(d) Growth factor g n under different T 


C. Numerical Result 

In the following, we show how to apply the above results 
to find the optimal n. We adopt WCDMA as the underlying 
cellular network technology. 100 omni-directional small cell 
base stations are deployed in a square domain of 1 km x 1 km. 
The network density is thus equal to: 

A = 10” 4 small base stations/m 2 . 

It is assumed that any cell synchronization identifier can 
be found in a radius Rb = 1 km. We take R m \ n equal to 2 
meters. The propagation path loss is modeled by the picocell 
path loss model l29l : 

PL [dB] (d) = 37 + 301og 10 (d) + 18.3/ ( 7TT-°- 46) , (23) 

where d is the distance from the base station in meters, / the 
number of penetrated floors in the propagation path. For indoor 
office environments, / = 4 is the default value (22]; however, 
here, the small cell network is assumed to be deployed in a 
general domain including outdoor urban areas where there are 
less penetrated walls and floors. So, we use / = 3 in our 
numerical study. 


It is assumed that the total transmission power including 
the antenna gain of each small cell base station is P Tx ,[dBml = 
32 dBm. Shadowing is modeled as a random variable with 
lognormal distribution with an underlying Gaussian distribu¬ 
tion of zero mean and 8 dB standard deviation. The signal 
strength received at any distance d from a base station i is 
expressible as: 

A,[dBm] (d) = Ar X , [dBm] ~ PL[dB](d) + . 

By (l23k we have: 


A,[dBm](rf) = ^Tx,[dBm] - 37 - 18 . 3 /^- 0 - 46 ) 

-301og 10 (d) + A; dB . (24) 


The parameters A and /3 appearing in l\ = A dr X, can 
be identified from (l24l) after converting the received signal 
strength from the dBm scale to the linear scale: 


The received noise power Nq is given by: 

Ao,[mW] = Ai^Kelvin X AAF[ W ] X W[Bz] x 10' 3 , 
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where the effective bandwidth Wy Hz ] = 3.84 x 10 6 Hz, 
is the Boltzmann constant, and Tkeivin is the temperature in 
Kelvin, fc B Tkeivi„ = 1.3804 x 1CT 23 x 290 W/Hz and NF [dB] 
is equal to 7 dB. 

It is assumed that the mobile is capable of scanning eight 
identified cells within 200 ms |28il . So, the average time 
needed to scan one cell is given by s = 25 ms. 

In order to check the accuracy of the approximations used in 
the analysis, a simulation was built with the above parameter 
setting. The interference field was generated according to a 
Poisson point process of intensity A in a region between R min 
and i?oo = 100 km. For a number n, the maximum of SINR 
received from n base stations which are randomly selected 
from the disk B between radii /f T11 jT1 and Rb was computed. 
After that the expectation of the maximal capacity E{£ n } 
received from the n selected BSs was evaluated. 

In Fig. 0 a), the expectation of the maximal throughput 
E {£„ } for different n is plotted, as obtained through the an¬ 
alytical model and simulation. The agreement between model 
and simulation is quite evident. As shown in Fig. 0a), E{£„} 
increases with n, though the increasing rate is slow down as 
n increases. Note that in Fig. [5j a), E{} is plotted after 
normalization by E{^ 25 o}. 

Fig. 0b) gives an example of acceleration p n for T = 0.5 
second and L(n) = n x 25 ms. In the plot, p„ is normalized 
by its maximum. Here, an agreement between model and 
simulation is also obtained. We see that p n first increases 
rapidly with n, attains its maximum at n = 42 by simulation 
and n = 43 by model, and then decays. 

Next, using the model we compute the optimal number of 
cells to be scanned and the growth factor g n for different T. 
Note that in ((Til . the factor T/(T + L(n )) can be re-written 
as: 

T = 1 [T< oo, 

T + L(n) l+nxs/T ° [L(n) = n x s 

It is clear that this factor also depends on the ratio T/s. 
Fig. 0c) plots the optimal n for different values of T/s. Larger 
T/s will drive the optimal n towards larger values. Since T 
can be roughly estimated as the mobile residence time in a 
cell, which is proportional to the cell diameter divided by the 
user speed, this can be rephrased by stating that the faster 
the mobile, the smaller T and thus the fewer cells the mobile 
should scan. 

Finally, Fig. 0d) plots the growth factor g n with different 
T. In Fig. 0d), the “limiting case” corresponds to the case 
when T —> oo or L(n) = 0. We see that g n is quite stable 
w.r.t. the variation of T. Besides, g n flattens out at about 30 
cells for a wide range of T. Therefore, in practice this value 
can be taken as a recommended number of cells to be scanned 
in the system. 

V. Concluding Remarks 

In this paper, we firstly develop asymptotic properties of the 
signal strength in cellular networks. We have shown that the 
signal strength received at the center of a ring shaped domain 
B from a base station located in B belongs to the maximum 


domain of attraction of a Gumbel distribution. Moreover, the 
maximum signal strength and the interference received from 
n cells in B are asymptotically independent as n —> oo. The 
above properties are proved under the assumption that sites are 
uniformly distributed in B and that shadowing is lognormal. 
Secondly, the distribution of the best signal quality is derived. 
These results are then used to optimize scanning in small cell 
networks. We determine the number of cells to be scanned for 
maximizing the mean user throughput within this setting. 
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Appendix A 
Proof of Lemma[2] 

Let di = y - x, be the distance from a site located at 
Xi £ K 2 to a position y £ R 2 . Under the assumption that site 
locations are uniformly distributed in B , the distance di from 
a site located in B , i.e., x^ £ B , to the center of B has the 
following distribution: 


F D (d) = P [dj < d\ 


nd 2 - 7 ri? 2 nin 

7tR 1 - ^ R Ln 


d 2 ~Rl in 


o2 _ c>2 

1X B ■H'min 


Let Ui = Ad i d, for f3 > 0 , its distribution is equal to: 

Fjj(u) = — c(u~^ — a - ^), for it £ [a, b), 


where c = A? (.- R 2 min ) \ a = AR/, and b = ARj[ 
The density of Ui is given by fu(u ) = (2 c/ f3)u~ l ~ 2 ^. 
Thus, the distribution F P of the power Pi is equal to: 


F P {x) = [ F x {-)fu{u)du. (25) 

J u—a ^ 

Substituting Fx with lognormal distribution of parameters 
(0 ,<Tx) and fu given above into (l25b . after changing the 
variable such that t = log(^), we have: 


F P {x) 



2 _ 

P 


log(f) 


e p di 


i°g(f) 


r log (t) 2 t . t 
e p erf (- 


'iog(f) 


k \/ 2 i 


ax 


)df) 


where the first integral is straightforward. By doing an integra¬ 
tion by parts of erf () and eP d t for the second integral, 
we get: 


Fp(x) 



2/3 

p — 

2 


e p 


— cl T \ 
+ e3erf (^—) 


— e p' 2 er: 


f( 


V%a x 


V?ax 
v /2 a x 

P 


log(f) 

*=l°g(§) 


After some elementary simplifications, we can obtain: 

-2/1 1 rA°& x ~ Mi ■ 


F F ( x ) = c{a-}(- + r vt(^±l)) 
(> + ‘ ert( >— 


- b~n 


V2, 


+e v x~^ 

1 


1 1 logx-/i 3 , 

* 2 ‘ 2 erf( ^T' 


1 loga;-^4 

5 T 5 erf( 


>]} 


where v = =^-, /o = log a, ^3 = Mi + 2a x /P, M 2 = log&, 
and fj ,4 = /j .2 + 2a\/p. Let Gj , j = 1,4, be the lognormal 
distribution of parameters (/ ij,ax ), j = 1,..., 4, Fp can be 
rewritten as 0 . □ 


Appendix B 
Proof of Lemma[3] 

Let Gj(x) = 1 — Gj{x) and note that c{a~~P — b~?) = 1, 
we have from 0 : 


Fp(x) = 1 — c^a PGi(x) — b pG 2 {x) 


2 - 


— e v x 0 G 3 (:e) + e":E ^G 4 (a;)j. 


This yields the tail distribution F p = 1 — Fp: 

F p (x) = c|a _ ^Gi(a:) — b~^G 2 (x) 


~?G 3 (:e) + e":r ^G 4 (a:)|. (26) 


For d26l >. we have Gj(x ) = \ erfc{(loga: — /Zj)/(v/2ox)}. 
An asymptotic expansion of erfc(x) for large x fl30l 7.1.23] 
gives us: 


G 3 (x) 


ax 


00 V / 27r(logx - /. 13 ) 
_ ax _ 

\/27r(log x — m - 


f (log X — /1 3 ) 2 'I 

exp l- 2% -J 

(log X — fl! - 


exp« 


2 a\ 


ax a e v x ? 

\/27r(log x - m) 


exp ] -- 


(lo g x /li ) 2 
2 a~ x 


1 


1 - 


2 ai 


0(iogx—m) 

in which after a Taylor expansion of the last term on the right- 
hand side, we can have: 


e v x p G 3 (a;) ~ a p Gi(x) + 

00 L 

This implies that 


2 a\ ex P ( - 


(log x—p,\) 2 




7T /3 (log X m) 2 


a PGi(x)—e v x ' 3 G 3 (x) 

/— 2 0 _ ( (logx-/xi ) 2 \ 

<2 a ?a\ ex P l- 2 pr - ) 

13 (log x fi \) 2 
In the same manner, we have 


b PG 2 {x) — e' y x pGa{x) 


(27) 


77,-2 o », rn / (Iog^~M2) 2 \ 

2b p a\ ex P l ) 

13 (log x fjj 2 ) 2 


(28) 
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A substitution of ( f271 ) and d28| ) into ( f26t results in 

F ( 2 .) ^ ^ ca x ex P ( ~ ~ ^) 2 /(2o-|)) 


7 T /3 l (log a: — M 2) 2 

_2 exp ( - (logx - Hi) 2 /(2cr 2 x )) 

a (logx- mi ) 2 


}• (29) 


Moreover, b > a yields M 2 — Mi = log( 6 /a) > 0- Then, we 
have the following result for large x: 

exp(- (log 2 y i) 2 )/(logx-Mi ) 2 


exp ( - (los 2 y 2)2 )/(logx - M 2) 2 


( log X M2 ^ 2 

Vlogx-uJ CXP l 2a 2 J 


_ 2 ,X "X -)• 0. 

2(7^- / x—foo 


- log x - Ml 

Taking this into account in ( |291 >. finally we have: 

_exp ( — (logx-M2) 2 /( 2cr D) 


F p (x) 


((logx - M 2 )/(v/2(t,y))^ 


2a/27t< 


naxn 


n = 


a x R min 


G 2 (x) _ 

logx -M2’ V^/3 R 2 B - R 2 min 


By ( l32l ). F q satisfies Lemma [ 6 ] with constants l = k, r = 
—2 ,77 = 1, and w = 2. So, Fq £ MDA(A) with the following 
normalizing constants: 

* (log n/rf)^~ l 1 1 

c n = -= Tt(logn) 2 , and 

LOT] 2 

J* _^logn^l/a; , t] 1 ^ r(log log n — log 77 ) + w log l ^ 
V / ' 2 - Nil 




(log n ) 1 i 


□ 


=( i„ gn) i + ILl2Sl25!i±l5eL. 

2 (log n) 2 

Then, by Lemma [5] we have 

lim P (M n < g(c^x + d*)) = A(x), xel. 

n—too \ / 

By a Taylor expansion of exp(v / 2<7xc*x), we have: 

lim P { e -(V2<7*d;+M) Mn < ! + y/ 2 a x c* n x + o«)} 

n—too K J 

= A(x). 

Since c* —>■ 0 when n —7 00 , we have 
M n - e V2<TX <+' i2 d 


V / 2CTA'C*e v/5crxd l* + ^ 2 


A, as n —> 00 . (34) 


Appendix C 
Proof of TheoremQ] 

We will use Lemma [3] and the following two lemmas to 
prove Theorem Q] 

Lemma 5 (Embrechts et al. m): Let Zi be i.i.d. random 
variables having distribution F, and T'n = max" =1 ZLet 
<7 be an increasing real function, denote Zi = g{Zf), and 
'I'n = max" =1 Z^ If F £ MDA(A) with normalizing constant 
c n and d n , then 

lim P(^ n < g(c n z + d n )) = A (z), z£l. 

n—f 00 

Lemma 6 (Takahashi l\3 If ): Let F be a distribution func¬ 
tion. Suppose that there exists constants u> > 0, l > 0, 77 > 0 
and rdl such that 


lim (l — F(x))/ (lx r e~ rlx “) = 1. 


(31) 


For m £ R and a > 0, let F* = F((x — p)/cr). Then, F* £ 
MDA(A) with normalizing constants c* = crc n and d* = 
ad n + p, where 


Cn 

dn 


(log n/ 77 ) 1 

w? 7 


and 


log 77 s i/ui 77 1 / aj r(log log n — log 77 ) -fwlog/ 

7? w 2 (logn) 1- - 


Let g(t) = e^ axt+fJ - 2 be a real function defined on R, g 
is increasing with t. Let Qi be the random variable such that 
Pi = 9 {Qi)- By ( 1 10ab of Lemma [3] the tail distribution Fq is 
given by: 

F q {x) = F P (e V2crxx+IJ - 2 ) ~ KX~ 2 e~ x \ as x -7 00 . (32) 


Substituting c* and d* from (f33l) into ( 1341 . we obtain c n and 
d n for O. The conditions f? max < 00 , > 0 and cr x > 

0 provide tt > 0. This leads to d n > 0, and consequently, 

c n > 0. □ 


Appendix D 
Proof of Lemma|4] 

Under the assumptions of the lemma, the interference field 
can be modeled as a shot noise defined on R 2 excluding the 
inner disk of radius F m in- Hence, using Proposition 2.2.4 in 
[|2Qj . the Laplace transform of I is given by: 

Ci(s) = exp | — 27tA f ^1 — Eje -3 ^ * }^rd?’|. (35) 

" ^min 

Noting that 


w) = Ci{—jw), wet, 

we have from (l35l > that: 

r ju.AXj 

4>i(w) = exp | — 27 tA / (1 — E{e ^ })rdr|. (36) 

“ -Rmin 

Using the change of variable t = \w\Ar~P, we obtain 


7*+OO 


jwAXj 

(1 — E{e r-P })rdr 


(A\w\) 2 /P 


•W 1 1 — E{e Jsign (“i tXi } 

min «- J 


i2//3+l 


df, (37) 


where E{e Jslsn ( u 'i tAi } = fx{sign(w)t). So, substituting this 
into (f36t . we get the first part of the Lemma 0] 
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From (fl5l >, for all p = 1,2,..., we have: 

\(t>i(w)\ p = exp ( 


A\w\ 

— p7rAa(A|w|) Q E{ f Rtnin - — cos (^-^t) ^ 

Jo 


t a+l 


Since 1 — cos {tXJ) > 0, Vi £ R, we have 

Afwl 

‘W~ 1 — cos (tXA 

min V L ' 


E 


t a+1 


di} > 0. 


Therefore 


\J>i(,w)\ p < exp(-c|iu|“), 


where 


h(v,u) 


(38) 


(39) 


1 if v e [u, Ytl u - 1 ] and u> 7 

0 otherwise 


It is easily seen that h(v, u) is square-integrable with respect 
to v, and its Fourier transform w.r.t. v is given by: 

- w 10 if zt < 7 

f J±i„-1 , (46) 

^ l e 3 dv if u > 7 


(40) which yields: 


where c is some positive constant, and hence the right hand- 
side of this is an absolutely integrable function. This proves 
the second assertion of Lemma 0] 

Under the assumption that AR~? n ks 00 , </>/ can be 
approximated by: 

4>i(w) ~ exp( — 7rAa(A|zn|)“ J -— 

(41) 

For 0 < a < 1, we have 

r°° l _ e jsign(w)tX, 

Jo 


t a+l 


df 


f*+oo 


1 - E{e- 


Jsign(w)tXi 


t a+1 


dt 




= -r(-a)|sign(zz;)A' i |“e- Jsign(u ' ) ^ (42) 

Since X t > 0, we can write |sign(zu)Aj| a = Xf. Taking 
expectations on both sides, we get 


.(^ e -j W u _ e j W (l-^u)^ ifti>7 - (47) 

Besides, according to Lemma 0] we have that (f>i £ L and 
4>i £ L 2 , where L 2 is the space of square integrable functions. 
And so, by Theorem 3 in l32l p.509], // is bounded continuous 
and square integrable. Hence, by applying the Plancherel- 
Parseval theorem to the inner integral of (l44l) . we have 

/•OO f‘OC 

/ h(v, u)fi(v)dv = / h(—w, u)fi(w)dw, (48) 

J 0 J —00 

where fi{w ) is the Fourier transform of //( v). Take (07} into 
account for (08]) and (1441 ). we have: 

Fv n ( 7)=/ {fM n {u)j h(—w, u)fi(w)dw\du, (49) 
J ■7 ^ J — OO ^ 


where we further have 


= -E{Xf}F(-a)e 

= cos(^)[l - jsign(w) tan®]. 

a 2 2 

Substituting this into (|4IT ) and noting that 

E{X,“} = exp (^a 2 a 2 x ), (43) 

for Xj lognormally distributed, we obtain ( fl 6 l> . □ 


h(—w, u)fi(w)dw 


P + OO 


2 n ./_„ * 2 jt 1 v 2 tt 


' —OO 
/*+ 


= 4 / 


2zr ,/g V 27T 27 t 2zr 2zr 


Note that 


w 




(51) 


Appendix E 
Proof of Theorem^] 

Under the assumption that sites are distributed as a homoge¬ 
neous Poisson point process of intensity A in B, the expected 
number of cells in B is tt\(R 2 b — i?X n ). We assume that 
n\(R B — f?min) * s muc h larger than n, which ensures that 
there are n cells in B with high probability, so that Y n is well 
defined. 

Under the conditions of Theorem 0] M n and I are asymp¬ 
totically independent according to Corollary 0] So, by substi¬ 
tuting (fJ~4t into ©, we have: 


P{Y n > 7 } 


And under the assumption that AR m ? n ss 00 , J>i is ap¬ 
proximated by (06}. Thus, by (06} and (07}, we have for 

w £ [ 0 , + 00 ): 


h( w ( w \ ~ e ^ f 
h (—’ u )fn-—) ~ ——v 


t j(-wu-\-dw a tan 


' 27T 


27T jw 


-j(—w+w u—Sw a tan 


} (52) 


and 


-5w° 




0 -j(—wu-\-5w a tan - 


fM n (u) h(v,u)fi(v)dvdu (44) 


27T 27T jw 

_|_ g R-W+W^-U-Sw 01 tan^)^ 
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By d52> and Q, we get 

1 \i, w . ; ,-w . - -w f w 

— & w 7TQ 

sin(—-wm + tan —) 


7TW 


+ sin (—w + W————U — Sw a tan 

7 2 ■ 

2e“ ,5,u '' 


■ sin in 


TTW 


t “ / \ 

(in—— 
v 2 7 ; 


X COS (wit + 


27 2 ' 


Substitute the above into (f50b and then into ( l49l >. we have 

O- □ 











